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Abstract: By exploiting a recently developed connection between Heun's differential equation and the 
generalized associated Lame equation, we not only recover the well known periodic solutions, but also obtain 
a large class of new, quasi-periodic solutions of Heun's equation. Each of the quasi-periodic solutions is 
doubly degenerate. 
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Heun's equation, a second order linear differential equation with four regular singular points has been 
extensively discussed in the mathematics literature [T][21[3|- In recent years, this equation has also appeared 
in a number of physical problems, like quasi-exactly solvable systems sphaleron stability [5], Calogero- 
Sutherland models jH], higher dimensional correlated systems jjj, Kerr-de Sitter black holes jSj, and finite 
lattice Bethe ansatz systems [D]. 

The so-called periodic (also termed as polynomial) solutions of Heun's equation have been well studied. 
However, as emphasized in ref. PQ, much less attention has been devoted to the quasi-periodic (also termed 
as non-polynomial) solutions. In this letter we obtain a large class of (mostly new) quasi-periodic solutions. 
We shall show that each such solution is doubly degenerate. 

The canonical form of Heun's equation is given by [Q 

G(x) = 0, (1) 



d 2 / 7 5 e \ d adx 
+ - + r + )^ + 



dx 2 \x x — 1 x — c/dx x(x — l)(x — c)_ 
where a, (3, 7, 5, e, q, c are parameters. The parameters are not all independent; the constraint relation is 

7 + «5 + e = a + /3 + l. (2) 

The four regular singular points of the differential equation are located at x = 0, 1, c (7^ 0,1), and the 
point at infinity. 

Let us make a change of independent variables using the transformation x = sn 2 (y,m) ^U], where 
m = 1/c. Then Heun's equation takes the form 



F"(y) + 



1 - 2c)m Sn ^' m ^ >Cn ^ ,m ' > + (1 - 25) sn (^' m ) dn (^' m ) + _ ^ cn (y , m )dn(y , m 



dn(y,m) cn(y,m) sn(y,m) J 

— [Amq — 4a/5msn 2 (y, m)]F(y) = , (3) 



where G{x) = F{y). The periodic solutions of eq. © correspond to the polynomial solutions of eq. Q 
while the quasi-periodic solutions of this equation correspond to the non-polynomial solutions of eq. 

Before describing the quasi-periodic solutions of eq. (j3J), it may be worthwhile explaining how we 
arrived at such solutions. A few years ago, we studied in some detail the Schrodinger equation for the 
associated Lame (AL) potentials I12j 

V(y, m) = a(a + l)msn 2 (y, m) + b(b + l)msn 2 (y + K(m),m) = a(a + l)msn 2 (y, m) + b(b + l)m — 77—^ — - , 

dn (y,m) 

(4) 



which after a transformation can be shown to be a special case of eq. ©• In particular, the band edges 
and the mid-band states of several of these potentials were studied, when a and b were related in various 
specific ways. Further, when a, b are both integers, we showed that these potentials had the special feature 
of possessing only a finite number of band gaps. This study has been extended to the generalized associated 
Lame (GAL) potentials [H Ej 

V(y, m) = a(a + l)msn 2 (y, m) + b(b + l)msn 2 (y + K(m),m) 

+ /(/ + l)msn 2 (y + K(m) + iK'(m),m) + g(g + l)msn 2 (y + iK'(m),m) 
= a(a + l)msn 2 (y, m) + b(b + l)m^4 + /(/ + i) dn ^ m ) + g(g + 1 



dn 2 (y,m) cn 2 (y,m) sn 2 (y,m) 

(5) 

The point to note is that after a transformation, the Schrodinger equation for the potential given in eq. 
© is identical in form to eq. (jSJ. In particular, let us start from the Schrodinger equation 

-^p^ + V(y,m)iP(y) = EiP(y), (6) 

with V(y,m) given by eq. ©. On substituting 

iP{y) = dn- b (y)cn~f(y) S n-3( y )<P(y) , (7) 

one can show that <t>{y) satisfies Heun's equation (j2J). More precisely, 4>{y) satisfies the differential equation 



0"(y)+2 



, sn(y, m)cn(y, m) „sn(y, m)dn(y, m) cn(y, m)dn(y, m) 
mb — r h/ -, r g- 



dn(y,m) ' cn(y,m) ' sn(y, m) 



0'{y)-[R-Qmsn 2 (y,m)](/)(y) = 0, 

(8) 

where 

R = -E + m{g + bf + (f + g)\ Q = (b + f + g)(b + f + g - 1) - a(a + 1) . (9) 

Thus, once we obtain solutions of the Schrodinger equation for the GAL potential ©, then we can imme- 
diately write the solutions of eq. © and hence that of the original Heun's eq. (0) with the identification 

Aa(3 = Q, 4mq = R, F(y) = <f>(y) . (10) 
A few comments are in order. 



1. The Schrodinger eq. © for the GAL potential © is invariant under y — > y + K(m) provided 
a and b are interchanged and so also are / and g. Hence, the eigenvalues E are invariant under 
a *-* 6 and / <-> g while the corresponding eigenfunctions are related to each other by the translation 
y —* y + K(m). Likewise, the eigenvalues are also invariant under a / and b ^ g, while the 
corresponding eigenfunctions are related to each other by the translation y —* y + K{m) + iK'{m). 
Similarly, the eigenvalues £7 are invariant under the transformation a <-> g and 6 <-> /, while the 
corresponding eigenfunctions are related to each other by the translation y — > y + iK'(m). 

2. The Schrodinger eq. © for the GAL potential © is also invariant under a — ► — o — 1 and/or 
6 — ► —6 — 1 and/or / — > — / — 1 and/or g — ► — g — 1. As a result the eigenvalues i£ as well as the 
corresponding eigenfunctions are invariant under one or several of these transformations. 

3. It may be noted that, except for the invariance under a — > —a — 1, eq. © is not invariant under 
any of the above transformations. However, the connection between the GAL problem and eq. © 
and the invariances of the GAL equation, can be exploited to obtain several more solutions of Heun's 
equation. For example, from eq. @, it follows that if under any of the above transformations, if 
bi, /i, <7i change to &2> f2, 92 and the energies E remain invariant, then the corresponding eigenvalues 
R = Amq of Heun's eq. (j3J) are related by 

Rl - m (b + 1 + 9l f - m{h +gi f = R 2 - m (b 2 + g 2 f - (f 2 + g 2 f . (11) 

Our strategy is now clear. We shall first obtain solutions of the Schrodinger eq. © for the GAL 
potential © and then using the connections given in eqs. H10|) and (Jllj) and the symmetries of the GAL 
equation, we shall obtain a host of solutions of Heun's equation. This strategy is demonstrated below by 
discussing one example in detail. As mentioned previously, the focus here is on the quasi-periodic solutions 
of the GAL equation. 

We now show that when either a + b + f + g and /or a — b — f — g is an arbitrary half- integer (but not an 
integer), then one can obtain doubly degenerate eigenstates of the GAL eq. © which correspond to the 
mid-band states (rather than band-edges) of this periodic problem. In the special case when a + b + / + g 
and/or a — b — f — g is an integer, then one obtains two distinct (non-degenerate) eigenstates which 



correspond to the band edge eigenstates. It turns out that depending on whether b or / or g is half- 
integral (while the other two are integral), we need to use a different ansatz. Let us consider all three cases 
one by one. 

Case 1: b half- integral 

We start from eq. © and substitute the ansatz 

cj)(y) = [cn(y, m) + isn(y, m)] t Z(y) , (12) 

where t is any real number. It follows that Z(y) satisfies the equation 

7"f \.ro-,M ( sn(y,m)cn(y,m) cn(j/,m)dn(j/,m) dn(y , m)sn(y , m) 

Z (y) + [2itdn(y, m) + 2mb — r 2g h 2/ \Z (y) 

an[y,m) sn{y,m) cn{y,m) 

, m , *2\ 2, \ n 'j. cn(y,m) sn(y,m) 



+ [-(R + r) + {Q + r)msxf(y, m) - 2itg ; ( + 2it/(l - m) 

sn(y, m) cn(y, m) 

+imt(2b + 2f + 2g- l)sn(y, m)cn(y, m)]Z(y) = , (13) 

where R and Q are as given by eq. ©• Not surprisingly, ^(y) = constant is a solution with energy 
E = (4t 2 + m)/4 provided / = g = 0, 6= 1/2, a = t - 1/2 (i.e. b + f + g = 1/2). 

One can build solutions for higher half-integer values of b + / + g from here. In particular, if b + / + g = 
2M + i, we choose the ansatz (M = 0, 1, 2, ...) 

M M-l 

Z (v) = X] Ak S1)2k (yi m ) + cn (2/' m )sn(y, m) .Bfe sn 2fc (y, m) , (14) 
fc=o fc=o 

while if b + / + g = 2M + 3/2, then we consider the ansatz (M = 0, 1, 2, ...) 

A/ M 

Z(y) = cn(y, m) ^ vl fc sn 2A: (y, m) + sn(y, m) B k sn 2k (y, m) . (15) 

A;=0 k=0 

Substitution into eq. (|T^|) followed by lengthy but straightforward algebra yields analytic expressions for 
the energy eigenvalues for arbitrary M and b = 1/2, 3/2. In particular, for b = 1/2, we find 

6 = 1/2, f + g = N, g=p, a = t-l/2, E = [t 2 + m(g + bf] , (16) 

where both /, g are nonnegative integers satisfying / + g = N with p,N = 0,1,2,... . Since p < N, 
it follows that for a given N, there are N + 1 different solutions. Similarly, when b = 3/2, we get 
9 =p,f = N -p,a = t- 1/2, and 



E = [1 + t 2 + m{g + b) 2 } - m(2g + 1) ± J(2c/ + l) 2 m 2 + 4m(iV + 1)(/ - g) + 4(1 - m)t 2 (17) 



where, / and g are again nonnegative integers. Note that, in this case too, for a given N, there are N + 1 
different solutions. In all cases, the corresponding eigenfunctions have the form given in eqs. (|14|) and 
(|15|) . For small values of N, the explicit coefficients A^^B^ appearing in the eigenfunction expressions can 
be easily written down. 

We can immediately write down the solutions of Heun's eq. © by making use of eq. (|1U|) . In particular, 
we find the following two classes of solutions: 



l = \-V^=\-N+p, e = 0, a = -^J, = -E+l, Amq = N 2 -t 2 ] (18) 

1 A 1 AT J- 1 N+l-t N+l+t 

7 = p , o = iv+p, e = — 1 , a = , p = , 

'2 2 2 2 



Amq = N — t — 1 + (2p + l)m ± ^ (2p + l) 2 m 2 + Am(N + l)(iV - 2p) + 4(1 - m)i 2 . (19) 

The corresponding eigenfunctions are generically as given by eqs. (|12[). I|14[) and l|15[) with F(y) = </>(y). 
As an illustration, for N = 1, the eigenfunction is Z(y) = Acn(y,m) + Bsn(y,m) with ^ = it for 
/ = i j5 = 0,6 = 1/2, while for g = 1, / = 0, b = 1/2, one gets f =i. 
Several remarks are in order at this stage. 

1. Since the Schrodinger eq. © for the GAL potential (JSJ as well as Heun's eq. (J3J) are invariant 
under y — > y + 2K{m)andy — > y + 2iK'(m), hence F(y), F(y + 2K(m)) and F(y + 2iK'(m)) 
are all eigenfunctions of Heun's eq. © with the same eigenvalue. As a consequence, F(y) = 
[cn(y, m) — isn(y, m)] t Z(y) is also an eigenfunction of Heun's eq. © with the same eigenvalue. Thus 
for any nonintegral t, each level is doubly degenerate. The same remark also applies to the other two 
solutions (when / or g are half integral) discussed below. 

2. For the special case of / = g = and t = 1/2, the results (as expected) are identical to those already 
obtained in ref. |12j . 

3. For integral t, both a, b are half integral (while /, g are integral) and each of these solutions reduces 
to two nondegenerate periodic solutions of the GAL equation. 
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4. We can generate new solutions of Heun's eq. © by considering the transformations y — > y + 
K(m),y —* y + iK'(m),y —* y + K(m) + iK'(m) and the fact that the eigenvalues are invariant 
provided we interchange a, 6, /, <? appropriately as explained above. Obviously, each of these new 
solutions is also doubly degenerate. In particular, by starting from solution and using eq. 
the three sets of new solutions of Heun's eq. © are 

1 P 1 N + t-2 

7 = - — N + p , o = - — p, e = 1 — t , a 



2 2 ' ' ' 2 

/3 = ^ — , 4mg = TV 2 - t 2 + m(N + t)[N + t - 2p - 1] , F = F(y + K(m)) . (20) 



1 N + t-p-1 
7 = l-t, 5 = 0, e = --iV + p,a = ^ , 

iV + i 

/9 = — , 4mg = m(N + t)[JV + t - 2p - 1] , F = F(y + iK'(m)) . (21) 



7 = 0, d = 1 - i, e = --p, a = , 

P = -^y^ , 4mg = , F = F{y + K(m) + ifir'(m)) . (22) 

Similarly, by starting from solution 1)19(1 with 6 = 3/2, three sets of new solutions can be immediately 
written down. 

5. For low values of N, the corresponding eigenfunctions can be explicitly shown. For example, for 
g = 0, / = 1, i.e. p = 0, N = 1, the eigenfunction corresponding to solution (JT8|) is given by 

F(y) = [cn(y, m) + isn(y, m)]* [cn(y, m) + itsn(y, m)\ , (23) 

while the eigenfunctions corresponding to solutions 1)20(1 to (|22|) are given by 

F = F{y + K{m)) oc [cn(y, m) + iy/l — msn(y, m)]*[tcn(y, m) + iy/l — msn(y, m)][dn(y, m)]~^ 1+ ^ , 

(24) 

F = F(y + iK'(m)) oc [1 - dn(y, m)]* [i - dn(y, m)\ [sn(y, m)]~ (1+t) , (25) 
F = F(y + K(m) +iK'(m)) oc [dn(y,m) - Vl - m]*[idn(y,m) - y 7 ! - m][cn(y, m)]~ (1+t) . (26) 
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6. We can generate even more new solutions of Heun's eq. © by considering the transformations 

b — > —b — 1, / — > —f — 1 and g — > —g — 1 either singly or in various combinations. As an illustration, 

by starting from solution (|18|) and using eq. we get the following seven sets of new solutions. 

1 . 1 N-t-2 
7 = 2 ° = 2 +P ~ e = 2 ' a = 2 ' 

(3 = - N + t ~\ 4m q = N 2 -t 2 -2m(2p-l), F = I® . (27) 
2 dn (y) 

1 A ^ N + 1 + 1 - 2p 
7 = p , = — h iV — p , e = U , a = , 

2 2 2 

/^ = - iV + 1 2 t " 2 ^ 4mg = iV 2 -t 2 -(2p-l)(2iV-2p+l), F = ^g^y ■ (28) 

3 . 1 A7 n 2p + t + l-iV 
7=2+^' = -+ ^- ^, e = 0, a = , 

P = - 2P + I ~ t - N ,*m q = N*-t>-(2 P + l)(2N-2 P -l), F = . (29) 

1 . 3 Ar n N + t + 3-2p N + 3-t-2p 
7 = 2~ p > 5 = - + N -p, e = 2, a = , /3 = , 

Amq = N 2 - t 2 - 2m(2p - 1) - (2p - l)(2iV — 2p + 1) , F = ^ 2(y) Ji^ 1(y) ■ (30) 

3 . 1 AT _ 2p + t + 3-iV 2p + 3-t-iV 
7=2+^' 5= 2 +?, ~ iV ' e = 2 ' a = 2^ ' ^ = 2 ' 

Amq = N 2 -t 2 + 2m(2p + 3) - (2p + 1)(2JV - 2p - 1) , F = 2/ ... ; : • (31) 



7 = — hp, o = — h A* — p, e = (J, a = , a = 

2 2 2 2 



4m, = A^ 2 — t 2 + A(N + 1) , F = m ^_^ m ^ 1(y) • (32) 

3 r 3 Ar n N + t + 4 o N + A-t 
7 =-+p, 5= - + N -p, e = 2, a = , 



2 2 2 2 

4mq = N 2 -t 2 + 2m(2p + 3) + A(N + l), F = - 2/ - , - =-TT7 - ■ (33) 



7. Starting from these seven solutions, we can generate further new solutions by considering the trans- 
formations y — > y + K(m), y — > y + iK'(m), y — > y + K(m) + iK'(m) and interchanging a,b,f,g 
appropriately. One can show that in this way, by starting from the solution (|18[) with 6 = 1/2 one 
has 20 independent solutions. 



8. There is one more remarkable symmetry associated with eqs. (|14[) and (|15|). Note that eq. (|15|) is 
invariant under i — > — i followed by i — > — z. Under this transformation, ansatz (|14l) becomes 

<p(y) = [cn(y, m) - isn(y, m)Y t Z(y) , (34) 

and hence it follows that solutions with the ansatz (|14j) and ()34j) as well as with the ansatz (p(y) = 
[cn(y,m) — isn(y, m)] t Z(y) are degenerate in energy (i.e. have the same value of R). This then gives 
us 12 additional solutions. 

9. Thus starting with the solution (|18|) . for b = 1/2, we can generate 32 new sets of quasi-periodic 
solutions of Heun's equation. In particular, 8 of these solutions have 7 = i — p, | + p, 8 = \ + 
p — N, I + N — p, e = 0,2; another 8 of these solutions have 7 = \ + p — N, | + iV — p, 5 = 
h — p, § + p, e = 1 ± t, another 8 have 7 = 0,2,<5 = l±i,e = ^— p, § + p and finally, another 8 have 
7 = 1 ± i, 5 = 0, 2, e = ^ +p — N, | + iV — p. The corresponding values of a, (3 are easily obtained by 
using a = a+ "' +s + t ~ 1 / 2 ; ^ = 7+<5+<; ~ a ~ 3//2 while 4m(7 can be computed by using eqs. (jTU|) and (fT%)) . 

10. Proceeding in the same way, by starting from the solution (j!9|) with b = 3/2, we also generate 32 
independent solutions in each of which 4mq takes two possible values. The corresponding values of 
7, 5, e are the same as in the 6 = 1/2 case except that the values 0, 2 are now replaced everywhere by 
— 1,3 respectively. 

Case 2: / half-integral 

We start from eq. © and substitute the ansatz 

<t>(y) = [dn(y, m) + iksn(y, m)] t Z(y) , (35) 
where t is any real number and k = y/m. It follows that Z(y) satisfies the equation 

Z"(y) + {2iktcn(y,m) + 2m6 ™fo™Mi/.™) _ 2 cnfonQdnfom) + dnfa, m)sn(y m) 

an(y,m) sn(y,m) cn(y, m) 



+[-(# + mr) + (Q - t 2 )msn 2 (y, m) - 2itkg — ^ — f - 2ifctfe(l - m) , 

sn(y,m) dn(y,m) 

+ilfet(2& + 2/ + 2 5 - l)sn(y, m)dn(y, m)]Z(y) = , (36) 
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where R and Q are as given by eq. ©. Not surprisingly, Z(y) = constant is a solution with energy 

E = {Amt 2 + l)/4 provided b = g = 0, /= 1/2, a = t- 1/2 (i.e. b + f + g = 1/2). 

One can build solutions for higher values of b + / + g from here. In particular, for b + f + g = 2M + 1/2, 

we consider the ansatz (M = 0, 1, 2, ...) 

M M-l 
Z (v) = H A k sn 2fc (y, m) + sn(y, m)dn(y, m) ^ S A sn 2fc (y, m) , (37) 
fc=0 fc=o 

while if b + f + g = 2M + 3/2 we take the ansatz (M = 0, 1, 2, ...) 

M M 
Z(y) = dn(y, m) ^ A fc sn 2fc (y, m) + sn(y, m) ^ B k sn 2k (y, m) . (38) 
fc=o fc=o 

Substitution into eq. (|3f))) yields analytic expressions for the energy eigenvalues and eigenfunctions for 
arbitrary M for / = 1/2 and / = 3/2. In particular, for / = 1/2, we find that 

f = l/2, g=p, b + g = N, a = t - 1/2 , E = [mt 2 + (g + ff] , (39) 

where both b,g are nonnegative integers satisfying b + g = N with JV = 0, 1, 2, ... . 
Similarly, when / = 3/2, a = t — 1/2, g = p,b + g = N we find that 



£l = [(1 + t 2 )m + (g + f) 2 ] - (2g + 1) ± ^ (2g + I) 2 + 4m(iV + 1)(6 - g) - 4m(l - m)£ 2 (40) 

where, b and 5 are again nonnegative integers. In all these cases, the corresponding eigenfunctions have 
the form as given above in eqs. 1)3 7|) and (|38[). For small values of N, the explicit coefficients ^4^, Bk in the 
eigenf unction expressions can be easily written down. 

Using eq. @, we can write down the corresponding solutions of Heun's eq. ©• They are given by 

1= l-p, 5 = 0, e = ^-N + p, a = -—Zl, / g = __±* ! Amq = m(N 2 - t 2 ) . (41) 

1 c -1 1 JV + l-t jV + l+J 
7=2~^' <5 = - 1 > e =2 _iV+ P' a = 2 ' ^ = 2 ' 



4mq = m{N 2 - if - 1) + (2p + 1) ± y (2p + l) 2 + 4m(iV + l)(iV - 2p) - 4m(l - m)t 2 . (42) 

Corresponding to each of these two solutions, we can again write down 32 independent sets of solutions as 
shown above in Case 1 when b was half-integral. 
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Case 3: g half- integral 

We start from eq. (JSJ and substitute the ansatz 

(j)(y) = [dn(y, m) + /ccn(y, m)] t Z(y) , (43) 
where t is any real number. It then follows that Z(y) satisfies the equation 

Z"(y) + [-2khm(y, m) + 2m6 ^%!^%H _ 2 °i(y,m)dn(y,m) + 2 dn(y,m)sn(y,m) 

dn(y,m) sn(y,m) cn(y,m) 

2/.. _a n^i, cn (y^ m ) n7.^ dn (y' m ) 



+[-12 + (Q + F)msn 2 (y, m) - 2ktb—^p — f - 2ktf- 

dn(y,m) cn(y,m) 

+fct(26 + 2f + 2g- l)cn(y, m)dn(y, m)]Z(y) = , (44) 

where 12 and Q are as given by eq. 0- Clearly, Z(y) = constant is a solution with energy E = (1 + m)/4 
provided 6 = / = 0, 5 = 1/2, o = t - 1/2 (i.e. 6 + / + g = 1/2). 

One can build solutions for higher values of b+f+g from here. In particular, when b+f+g = 2M + 1/2, 
we consider the ansatz (M = 0, 1, 2, ...) 

M M-l 
Z (y) = X! A ksn 2k (y, m) + cn(y, m)dn(y, m) ^ B k sn 2k (y, m) , (45) 
fc=o fc=o 

while if b + / + g = 2M + 3/2 then we consider the ansatz (M = 0, 1, 2, ...) 

M M 

Z(y) = cn(y, m) A k sn 2k (y, m) + dn(y, m) B k sn 2k (y, m) . (46) 
fc=o fc=o 

Substitution into eq. Q44|) leads to analytic expressions for the energy eigenvalues and eigenfunctions for 
arbitrary M when b = 1/2, 3/2. In particular, for b = 1/2, we find that 

5 = 1/2, b + f = N, o = t-l/2, £= [(/ + 5 ) 2 + m( 5 + 6) 2 ], (47) 

where both 6,/ are nonnegative integers satisfying b + f = N with iV = 0,1,2,.... Similarly, when 
5 = 3/2, a = t - 1/2,6 + / = N we obtain 



£ = (/ + g) 2 + m( 5 + 6) 2 - [1 + 2f + (26 + l)m] ± ^/(l - m)[(2/ + l) 2 - (26 + l) 2 m] + 4mt 2 . (48) 
In all these cases, the corresponding eigenfunctions have the form given above in eqs. (|43|) and l|4fij). 
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We now write down the solutions of Heun's eq. (J3J) corresponding to solutions (|47j) and (|48|): 

7 = 0, 5 = ^-7V + p, e=^-p, a = -^ = ^, /3 = _ ^^' = ° ' (49) 

i r ^ /v _i_ 1 N + l-t N+l+t 

7 = — 1 , o = iV+p, e = p , a = , p = , 

' 2 ^'2 2 2 



4mg = 2(iV - p) + (2p + l)m ± \J {I - m)[(2N -2p + l) 2 - (2p + l) 2 m] + 4mi 2 . (50) 

Corresponding to each of these two solutions, we can again write down 32 independent sets of solutions as 
done above when b was half-integral. 

Summarizing, in this letter we have obtained new 192 sets of quasi-periodic solutions of Heun's eq. 
each of which is doubly degenerate. In each set, solutions exist with polynomials of arbitrary order N, and 
for each N, there are N + 1 distinct solutions. 

We would like to thank K. Takemura for informing us about some related work and recent developments 
in the mathematics literature on GAL potentials and Heun's equation. 
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